package com.fjzx.androidframe.global.util;

/**
 * 地图坐标器转换工具
 *
 * @author xiongyongcong
 */
public class XCoordUtil {
    /**
     * 圆周率PI
     */
    private static double PI = Math.PI;
    /**
     * 卫星椭球坐标投影到平面地图坐标系的投影因子
     */
    private static double AXIS = 6378245.0;
    /**
     * 椭球的偏心率(a^2 - b^2) / a^2
     */
    private static double OFFSET = 0.00669342162296594323;
    /**
     * 圆周率转换量
     */
    private static double X_PI = PI * 3000.0 / 180.0;

    /**
     * 百度地图坐标系(bd09)转国家测绘局坐标系(gcj02)
     *
     * @param latitude   纬度
     * @param longtitude 经度
     * @return double[] 0:latitude(纬度), 1:longtitude(经度)
     */
    public static double[] BD09ToGCJ02(double latitude, double longtitude) {
        double x = longtitude - 0.0065;
        double y = latitude - 0.006;
        double[] latlng = new double[2];
        double z = Math.sqrt(x * x + y * y) - 0.00002 * Math.sin(y * X_PI);
        double theta = Math.atan2(y, x) - 0.000003 * Math.cos(x * X_PI);
        latlng[0] = z * Math.sin(theta);
        latlng[1] = z * Math.cos(theta);
        return latlng;
    }

    /**
     * 国家测绘局坐标系(gcj02) 转 百度地图坐标系(bd09)
     *
     * @param latitude   纬度
     * @param longtitude 经度
     * @return double[] 0:latitude(纬度), 1:longtitude(经度)
     */
    public static double[] GCJ02ToBD09(double latitude, double longtitude) {
        double x = longtitude;
        double y = latitude;
        double[] latlng = new double[2];
        double z = Math.sqrt(x * x + y * y) + 0.00002 * Math.sin(y * X_PI);
        double theta = Math.atan2(y, x) + 0.000003 * Math.cos(x * X_PI);
        latlng[0] = z * Math.sin(theta) + 0.006;
        latlng[1] = z * Math.cos(theta) + 0.0065;
        return latlng;
    }

    /**
     * 百度地图坐标系(bd09) 转 地球坐标系(wgs84)
     *
     * @param latitude   纬度
     * @param longtitude 经度
     * @return double[] 0:latitude(纬度), 1:longtitude(经度)
     */
    public static double[] BD09ToWGS84(double latitude, double longtitude) {
        double[] latlng = BD09ToGCJ02(latitude, longtitude);
        return GCJ02ToWGS84(latlng[0], latlng[1]);
    }


    /**
     * 地球坐标系(wgs84)  转  百度地图坐标系(bd09)
     *
     * @param latitude   纬度
     * @param longtitude 经度
     * @return double[] 0:latitude(纬度), 1:longtitude(经度)
     */
    public static double[] WGS84ToBD09(double latitude, double longtitude) {
        double[] latlon = WGS84ToGCJ02(latitude, longtitude);
        return GCJ02ToBD09(latlon[0], latlon[1]);
    }

    /**
     * 地球坐标系(wgs84) 转 国家测绘局坐标系(gcj02)
     *
     * @param latitude   纬度
     * @param longtitude 经度
     * @return double[] 0:latitude(纬度), 1:longtitude(经度)
     */
    public static double[] WGS84ToGCJ02(double latitude, double longtitude) {
        double[] latlon = new double[2];
        if (outOfChina(latitude, longtitude)) {
            latlon[0] = latitude;
            latlon[1] = longtitude;
            return latlon;
        }
        double[] deltaD = transform(latitude, longtitude);
        latlon[0] = latitude + deltaD[0];
        latlon[1] = longtitude + deltaD[1];
        return latlon;
    }

    /**
     * 国家测绘局坐标系(gcj02) 转 地球坐标系(wgs84)
     *
     * @param latitude   纬度
     * @param longtitude 经度
     * @return double[] 0:latitude(纬度), 1:longtitude(经度)
     */
    public static double[] GCJ02ToWGS84(double latitude, double longtitude) {
        double[] latlon = new double[2];
        if (outOfChina(latitude, longtitude)) {
            latlon[0] = latitude;
            latlon[1] = longtitude;
            return latlon;
        }
        double[] deltaD = transform(latitude, longtitude);
        latlon[0] = latitude - deltaD[0];
        latlon[1] = longtitude - deltaD[1];
        return latlon;
    }


    public static double[] transform(double latitude, double longtitude) {
        double[] latlng = new double[2];
        double dLat = transformLat(longtitude - 105.0, latitude - 35.0);
        double dLon = transformLng(longtitude - 105.0, latitude - 35.0);
        double radLat = latitude / 180.0 * PI;
        double magic = Math.sin(radLat);
        magic = 1 - OFFSET * magic * magic;
        double sqrtMagic = Math.sqrt(magic);
        dLat = (dLat * 180.0) / ((AXIS * (1 - OFFSET)) / (magic * sqrtMagic) * PI);
        dLon = (dLon * 180.0) / (AXIS / sqrtMagic * Math.cos(radLat) * PI);
        latlng[0] = dLat;
        latlng[1] = dLon;
        return latlng;
    }

    public static double transformLat(double x, double y) {
        double ret = -100.0 + 2.0 * x + 3.0 * y + 0.2 * y * y + 0.1 * x * y + 0.2 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(y * PI) + 40.0 * Math.sin(y / 3.0 * PI)) * 2.0 / 3.0;
        ret += (160.0 * Math.sin(y / 12.0 * PI) + 320 * Math.sin(y * PI / 30.0)) * 2.0 / 3.0;
        return ret;
    }

    public static double transformLng(double x, double y) {
        double ret = 300.0 + x + 2.0 * y + 0.1 * x * x + 0.1 * x * y + 0.1 * Math.sqrt(Math.abs(x));
        ret += (20.0 * Math.sin(6.0 * x * PI) + 20.0 * Math.sin(2.0 * x * PI)) * 2.0 / 3.0;
        ret += (20.0 * Math.sin(x * PI) + 40.0 * Math.sin(x / 3.0 * PI)) * 2.0 / 3.0;
        ret += (150.0 * Math.sin(x / 12.0 * PI) + 300.0 * Math.sin(x / 30.0 * PI)) * 2.0 / 3.0;
        return ret;
    }

    public static boolean outOfChina(double latitude, double longtitude) {
        if (longtitude < 72.004 || longtitude > 137.8347) {
            return true;
        }
        if (latitude < 0.8293 || latitude > 55.8271) {
            return true;
        }
        return false;
    }
}
